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Abstract 

Based on previous studies of hard exclusive leptoproduction of pions 
in which the essential role of the pion pole and the transversity gener¬ 
alized parton distributions (GPDs) has been pointed out, we present 
predictions for the four partial cross sections of the exclusive Drell-Yan 
process, 7r~p —>• l~Nn. 

1 Introduction 

In recent years hard exclusive leptoproduction of mesons and photons have 
been studied intensively by both experimentalists and theoreticians. It be¬ 
came evident in the course of time that within the handbag approach which 
is based on QCD factorization in the generalized Bjorken regime of large 

^Email: goloskkv@theor.jinr.ru 
^Email: kroll@physik.uni-wuppertal.de 



photon virtuality and large photon-proton center-of-mass energy but hxed 
x-Bjorken, it is possible to interpret these processes in terms of generalized 
parton distributions and hard perturbatively calculable subprocesses with, 
however, occasionally strong power corrections for meson production (for a 
recent review see [1]). Exploiting the universality property of the GPDs, 
one may use the set of GPDs extracted from meson leptoproduction, in the 
calculation of other hard exclusive processes. Of particular interest are pro¬ 
cesses with time-like virtual photons. Thus in [2] predictions for time-like 
DVGS {'jp —)■ l~l^p) have been given, their experimental examination is still 
pending. The high-energy pion beam at J-PARG put into operation in the 
near future, offers the possibility of measuring another exclusive process with 
time-like virtual photons, namely the exclusive limit of the Drell-Yan pro¬ 
cess, n~p —)■ The purpose of this letter is to present predictions for the 

cross sections of this process taking into account what has been learned in 
the analyses of pion leptoproduction [3, 4]. The data on the cross section for 
leptoproduction [5, 6] demonstrate the prominent role of the contribution 
from the pion pole at small invariant momentum transfer, t, and it became 
evident that it is to be calculated as an one-particle-exchange (OPE) term 
rather than from the GPD E [7]. In the latter case the pion-pole contri¬ 
bution to the 7 r~^ cross section is underestimated by order of magnitude. A 
second important observation has been made in [3, 4]: The interpretation of 
the transverse target spin asymmetries in vr’*' leptoproduction measured by 
the HERMES collaboration [8] necessitates contributions from transversely 
polarized photons which are to be modelled by transversity GPDs within 
the handbag approach. This observation is supported by a recent GLAS 
measurement of tt® leptoproduction [9]. 

Since for the process 7 r~p —)■ the same GPDs contribute as for pion 

leptoproduction and the corresponding subprocesses are just s -H- it crossed 
ones ^ 

( 1 ) 

where s and u denote the subprocess Mandelstam variables, one can exploit 
the knowledge acquired there. One thus gains predictive power, there is 
no free parameter or soft hadronic matrix element left for the Drell-Yan 
process. Our analysis markedly differs from a previous study performed by 
Berger et al [11] where only predictions for the longitudinal cross section at 

detailed discussion of the space- and time-like connection of the leading-twist am¬ 
plitudes can be found in [10] 
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leading-twist accuracy has been given. It should be stressed that their and 
our predictions for that cross section differ by about a factor of 40 due to 
the different treatment of the pion pole contribution. Our hndings may be 
of help in the preparation of an Drell-Yan experiment [12]. Future data on 
the exclusive pion-induced Drell-Yan process may reveal whether or not our 
present understanding of hard exclusive processes in terms of convolutions 
of GPDs and hard subprocesses also holds for time-like photons . This is a 
non-trivial issue because the physics in the time-like region is complicated 
and often not understood. Thus, for instance, there is no explanation of 
the time-like electromagnetic form factors of hadrons [13]. Even the semi- 
inclusive Drell-Yan process was difficult to understand. It took a long time 
before the discrepancy between the theoretical predictions and experiment, 
known as the i^-factor, has been explained as threshold logarithms [14, 15] 
representing gluon radiation resummed to next-leading-log (NLL) accuracy. 


2 The handbag approach 


Here, in this section, we recapitulate the handbag approach. For more details 
of it we refer to our previous work [3, 4]. The process n~p —)■ is depicted 

in Fig. 1. We work in a center-of-mass frame in which p-fp' points along the 
positive 3-axis and we consider the kinematical range of large Mandelstam s 
(= (p -|- qY) and large photon virtuality ^ , Q'^, but small 


s — m 


( 2 ) 


the time-like analogue of Bjorken-x {m being the mass of the nucleon). Hence, 
skewness, dehned as 


p+ — p'+ r 
p+ -f p'+ 2 — r ’ 


( 3 ) 


is also small. 

Assuming factorization we can express the helicity amplitudes for 7 r~p 


"‘The Q'^-regions of quarkonia states have to be excluded. 
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Figure 1: The exclusive Drell-Yan process. The symbols in brackets denote 
the momenta of the respective particles. 


7 *n in terms of convolutions of GPDs and hard subprocess amplitudes 
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Explicit helicities are labeled by their signs or by zero, Cq denotes the positron 
charge and t' = t — Iq where to = —is the minimal value t 
corresponding to forward scattering. Terms of order t/Q'"^ are neglected 
throughout. The amplitudes for negative helicity of the initial state proton 
are obtained from the set of amplitudes (4) by parity conservation. The 
residue of the pion pole is given by 

Ptt = V^5'7rAf7vF7rArAr(t)Q^^FV(Q'^) (5) 

where Qt^nn (= 13.1 ± 0.3) is the familiar pion-nucleon coupling constant 
and En-ATAf is a form factor that describes the t-dependence of the coupling 
of the virtual pion to the nucleon. The pion mass, ttItt, is neglected except 
in the pion propagator. As we mentioned in the introduction we treat the 
pion pole as an OPE term. Therefore the full time-like electromagnetic form 
factor occurs in (5). Calculating the pion pole contribution from the GPD 
E as it is done in [11], one obtains to same the expression for it but with the 
leading-order (LO) perturbative result for the pion form factor. In (4) it is 
also allowed for a possible non-pole (n.p.) part of E. 
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For incident tt” mesons the p ^ n transition GPDs are required which, as 
a consequence of isospin invariance, are given by the isovector combination 
of proton GPDs [7] 

;y(3) = (6) 

The convolutions of the GPDs and the amplitudes "H for the subprocess 
Ti~q —)■ 7 *g read [3, 4] 

(77(3)) ^ j . (7) 


The helicity of the hnal state quark is A = /i + 1/2 with the photon helic- 
ity, p., being either zero or -1. Thus, the asymptotically leading longitudinal 
amplitude is related to a helicity-non-flip subprocess amplitude while, for 
transverse photons, a helicity-flip amplitude is convoluted with the transver- 
sity GPDs Ht and the combination Et = 2Ht + Et- As made explicit in 
(4) the transverse amplitudes are suppressed by p,T,/Q' as compared to the 
longitudinal ones. The mass parameter is related to the chiral condensate 


Pn 

rriu + md 


( 8 ) 


{rriui md are current quark masses). The subprocess amplitudes are calculated 
to LO of perturbation theory retaining quark transverse momenta, k^, and 
taking into account Sudakov suppressions while the emission and reabsorp¬ 
tion of partons by the nucleon happens collinearly to the nucleon momenta. 
This so-called modihed perturbative approach turns into the leading-twist 
result [11] for Q''^ —)■ oo. 

Since the Sudakov factor, exp [—S'], comprises gluonic radiation, resummed 
to all orders of perturbation theory in NLL approximation [16] which can only 
be efficiently performed in the impact parameter space, canonically conju¬ 
gated to the k^-space, one is forced to work in the b-space. Hence, 


= J dz(fb'^_x+{z,-h)F^x,o+{x,^,z,Q'‘^,h) 

X Q;s(/iij)exp [-S(x,b, Q'^)] . (9) 

The Fourier transforms of the hard scattering kernel and the light-cone wave 
function of the pion are denoted by F and T, respectively. The momentum 
fraction of the helicity -1-1/2 quark entering the pion is denoted by x; the 
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helicity of the antiquark is —A. For the renormalization scale we choose 
Hr = max(2;(5', (1 — z)Q', 1/6) and the factorization scale is 1/6. Following 
Li and Sterman [16] we only retain the most important quark transverse 
momenta which appear in the denominators of the parton propagators in the 
hard scattering kernels. Therefore, we can use the light-cone projector of a 
qq pair on an ingoing pion in collinear approximation [17] 


V 

• TT 


/tt 75 




/^7I 


$ 




(Z fjtu 


q^rf 
q ■ n 


^a{.z)-q^^ 


d^a{z) 

dz 



( 10 ) 


and replace the distribution amplitudes by light-cone wave functions. In (10) 
fn{= 132 MeV) is the pion decay constant, Nc the number of colors and n is 
a light-like vector which in a frame where the massless pion moves along the 
^(-direction is n = [0,1, 0j_]. Three-particle conhgurations, qqg, are neglected. 
Dirac, flavor and color labels are omitted for convenience. The hrst term in 
(10) is the well-known twist-2 part which is employed in the calculation of 
^o±,o+- Foi’ the accompanying light-cone wave function we take 


^ exp [-a^ki/(;^(l-;.))] (11) 

Jtt 

with the transverse size parameter = [v^tt/tt] ^ hxed from —?■ 77 

decay [18]. The twist-3 part of (10) is utilized in the calculation of Ti _ 

As the calculation reveals this subprocess amplitude is dominated by the 
contribution from <hp while the tensor term provides a correction of order 
t/Q'"^ which is neglected for consistency. For the wave function associated to 
$P (= 1), we use ^ 






/iraplkpl exp [—Opk^ . 


( 12 ) 


For the transverse size parameter, ap, we take 1.8 GeV 

^ Since quark and antiquark forming the pion have the same helicity, it may seem 
appropriate to use a h = ±1 wave function (for a particle moving along the z-direction). 
Such a wave function has been proposed in [19]. It is proportional to = fcj] ± ik^. Its 
collinear reduction leads to the tensor piece in (10), the important term ~ <I>p is lacking. 
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3 Predictions for the partial cross sections 


Before we present our predictions for the exclusive Drell-Yan process we 
specify the various parameters and soft hadronic functions we use in the 
evaluation. The form factor appearing in (5), is parametrized as 

_ A%-ml 



with An = 0.44 ± 0.04 GeV. For the time-like pion electromagnetic form 
factor also occuring in (5), we take the average of the data from CLEO [20] 
and BaBar [21] as well as a value derived from the J/T ^ 7r’''7r“ decay [22] 


Q'2|F^(g'2)| = 0.88 ± 0.04 GeV^ . (14) 


For its phase, exp [iS^Q'"^)], we rely on a recent dispersion analysis [23] which, 
for 2 GeV^ GeV^, provides 


5 = I.OMtt + 0.195(g'VGeV2 -2) - 0.029(g'VGeV2 -2)^ . (15) 


In the absence of any other information on this phase we use this parametriza- 
tion up to ~ 8.9 GeV^ where 5 = tt. For larger values of Q'^ we take 6 = n, 
the asymptotic phase of the time-like pion form factor obtained by analytic 
continuation of the perturbative result for the space-like form factor [13]. 

The GPDs are constructed with the help of the familiar double distribu¬ 
tion ansatz from the zero-skewness GPDs which are parametrized as ® 

K{x,^ = 0,t) = k{x) exp [t{b + a'In x] (16) 

where the forward limit, k{x), is an appropriate parton distribution (PDF) 
or is parametrized like a PDF with parameters htted to experiment. The 
GPD H{x,^ = 0,t) (including an error estimate) is taken from the recent 
analysis of the nucleon form factors [24] which, for this GPD, is based on the 
DSSV polarized PDFs [25]. A non-pole contribution to E is neglected, there 
is no clear signal for it in the data on pion leptoproduction. For the zero¬ 
skewness transversity GPDs, Ht and Et, the actual values of the parameters 
are specihed in [26]. They are oriented on lattice QGD results [27, 28] and 
lead to fair fits of the pion leptoproduction data [5, 8, 9] as well as of the spin 

®A more complicated profile function is adopted for H, see [24]. 
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Figure 2: Definition of the angles 0 and 6 . The latter angle is dehned in the 
rest frame of the virtual photon. 


density matrix elements and transverse target spin asymmetries for vector 
mesons [29, 26]. The errors of the transversity GPDs are estimated from the 
p-pole hts presented in [27, 28]. For the mass parameter (8) that controls 
the strength of the twist-3 amplitudes, we adopt the value ^ /ijr = 2 GeV 
valid at the scale 2 GeV. For its error we choose -1-0.55 and —0.15 [22]. The 
QGD coupling constant, Os, is evaluated from the one-loop epression for four 
flavors and Aqcd = 182 MeV. The time-like Sudakov factor is unknown, the 
continuation from the space-like to the time-like region is not well understood 
(see [30]). The replacement of by —Q'^ (see [30, 31]) leads to an oscillating 
phase but it is unclear whether these oscillations are physical or not. We 
therefore follow Gousset and Fire [30] and use the space-like Sudakov factor, 
as utilized in our previous work, also in the time-like region (with —)■ 

As shown in [32], for Q'^ less than 10 GeV^ the Sudakov fator is always close 
to unity except near b = 1/Aqcd where it drops to zero sharply. With the 
exception of this region the wave function provides the main suppression. 
Therefore, the detailed behavior of the Sudakov factor is not very important. 

The four-fold differential cross section for n~p —)• l~l^n read 




sin (26*) cos 0 daiT 


y/2 dtdQ’'^ 


where the angles 0 and 6*, specifying the directions of the leptons, are dehned 
in Fig. 2. The partial cross sections are related to the 7r“p —)■ 7 *n helicity 

’^According to the recent particle data tables [22] is rather 2.6 GeV. Using this 
value the normalizations of Ht and Et have to be altered accordingly since the fit to the 
pion leptoproduction data fixes the product of p,,r and the transversity GPDs. 













amplitudes (4) by 

ddL 

dtdQ''^ 

d(jT 

dtdQ'"^ 

da LT 
dtdQ'"^ 


^ |A^oi/',o+P, 

I/' 

/j,=±i,i^' 

K R-e^^ 0+(-^+i^',0+ ~ M.-v'fl+)\ ; 

v' 


daxT 

dtdQ'"^ 

The normalization factor 


Re^ 

v' 

reads (lepton masses are neglected) 
_ Clem 1 

dSTT^ (s — m^Y‘Q'‘^ 


(18) 


(19) 


In Fig. 3 we show our predictions for daL/dtdQ'"^ at Q'^ = 4 GeV^ and 
s = 20 GeV^ and daL/dQ'‘^ integrated over t' from 0 to —0.5 GeV^. The 
longitudinal cross section is heavily dominated by the contribution from the 
pion pole, that one from H, including its interference with the pion pole, 
amounts only to about 10% in the kinematical range of interest. The full 
result is markedly larger than our leading-twist result which is of the same 
order as that one quoted in [11]. This amplification is due to the use of 
the experimental value of the pion form factor (14) instead of its leading- 
twist result (~ 0.15 GeV^). We stress that the OPE contribution from the 
pion pole does neither rely on QGD factorization nor on a hard scattering. 
It is therefore not subject to evolution and higher-order perturbative QGD 
corrections. Because of the dominant contribution from the pion-pole and 
since we only consider a small range of Q'^ around 4 GeV^ the evolution of 
the GPDs is insignihcant and is therefore neglected. As opposed to [11] our 
interference term is positive. It is generated by the imaginary parts of (FT) 
and the pion-pole contribution while, in a LO leading-twist calculation, it is 
evidently under control of the corresponding real parts. Gonstructing H from 
the polarized PDFs derived in [33] instead from the DSSV ones [25] alters 
the predictions for the longitudinal cross section by less than the estimated 
errors displayed in Fig. 3. 

The transverse cross section is shown in Fig. 4. It is substantially smaller 
than the longitudinal cross section but much larger than the leading-twist 
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Figure 3: The longitudinal cross sections daL/dtdQ''^ (left) at Q'^ = 4 GeV^ 
versus t' and dai/dQ'"^ (right) versus Q'^. The thin solid lines with error 
bands represent our full results at s = 20 GeV^, the thick dashed ones those 
at 30 GeV^. The thick solid (dotted, thin dashed) line is the interference 
term (contribution from leading twist). The latter two results are 

multiplied by 10 for the ease of legibility. 


result. The uncertainty of our predictions is rather large and asymmetric 
due to the asymmetric error of The transverse cross section can be 
decomposed as (cf. (18) and (4)) ® 

= K[|V(-,o+t + 2|A1++,„+(^)r + . (20) 

The hrst term in (20), being related to the GPD Ht, is displayed in Fig. 4 
separately; it dominates this cross section. The second term, the pion-pole 
contribution, is rather small; it generates the little difference between the 
contribtuion from Ht and the full result for dax- The contribution from Et 
is tiny. 

The longitudinal-transverse interference cross section is also shown in Fig. 
4. The width of its error band is about a half of that of the transverse cross 

®The Et (tt) term behaves as behaves as a natural (unnatural) parity exchange while 
the Ht has no specific parity behavior [3, 4]. 
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Figure 4: The transverse cross sections daT/dtdQ''^ (left) at Q'^ = 4 GeV^ 
versus t' and dcTT/dQ''^ (right) versus Q'^. The thin solid dashed lines with 
error bands represent the full result at s = 20 GeV^, the thick dashed ones 
those at 30 GeV^ while the dotted line is the contribution from Ht- The thick 
solid line represents the longitudinal-transverse interference cross section. 


section, dair can be written as 
daiT 


dtdQ' 


1/2 


— K Re [2A4 q^ Q^A^++o+(vr) A4g_ g_|_A^—^o+] • 


( 21 ) 


Both the terms signihcantly contribute to dair- The transverse-transverse 
interference cross section is given by 


. ( 22 ) 

This cross section is very small. For instance, at Q'^ = 4 GeV^ and s = 
20 GeV^ it is less than 0.3 pb/ GeV"^. 

The cross sections decrease with growing s. As an example we show re¬ 
sults at s = 30 GeV^ in the plots. At, say, s ~ 360 GeV^ as is available from 
the pion beam at GERN, the longitudinal cross section is about 30 fb/ GeV^ 
at Q'^ = 4 GeV^. This is likely too small to be measured. 
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4 Conclusions 


We calculated the partial cross sections for the exclusive Drell-Yan process, 
'K~p —)■ within the handbag approach. In contrast to a previous 

study of this process [11] we treat the pion pole as an OPE term and take 
into account transversity GPDs. The parametrizations of the GPDs H, Ht 
and Erp as well as the values of other parameters appearing in the present 
calculation are taken from previous work [3, 4, 24], The generalization of our 
approach to K~p —)■ is straightforward. 

Future data on n~p —)■ l~l^n measured at J-PARG may allow for a test of 
factorization of the process amplitudes in hard subprocesses and soft GPDs. 
In contrast to pion leptoproduction where there is a rigorous proof for factor¬ 
ization of the amplitudes for longitudinally polarized photons, factorization of 
the exclusive Drell-Yan process is an assumption although it seems plausible 
that the factorization arguments also hold for time-like photons. However, 
Qiu [34] conjectured that factorization may be broken for the exclusive Drell- 
Yan process. If however factorization holds to a sufficient degree of accuracy 
future data on the exclusive Drell-Yan process may improve our knowledge 
of the GPDs. 

The exclusive Drell-Yan process also offers the opportunity to check the 
dependence of the TTTry vertex on the pion virtuality by comparing data on 
the time-like form factor measured in /+/“ —)■ with parametrizations 

of Tr^TT’*'* as part of the Drell-Yan analysis. The extraction of the 

space-like form factor from Ip data may beneht from that check. 
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